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ABSTRACT 

Over the last few decades, observational surveys have revealed that high-order 
multiple-star systems (e.g. triples, quadruples, etc.), and triples in particular, are com- 
mon in our Galaxy. In this paper, we consider the dynamical significance of this trans- 
formation in our understanding of stellar multiplicity. Using empirically constrained 
binary and triple fractions in those star clusters for which these values are available 
in the literature, we compare analytic rates for encounters involving single, binary, 
and triple stars. Our results show that, even for relatively low triple fractions, dy- 
namical interactions involving triples occur roughly as often as encounters involving 
either single or binary stars alone, particularly in low-mass star clusters. More gen- 
erally, using empirically-derived multiple star catalogues for the young star-forming 
association Taurus- Auriga and the Galactic field, we show that the data are consistent 
with the gravitationally-focused cross section for encounters increasing with increas- 
ing multiplicity. Consequently, triple stars, and even higher-order multiples, could be 
more important than previously realized for a number of astrophysical phenomena, in- 
cluding the formation and destruction of compact binaries and various types of stellar 
exotica, and the dynamical evolution of star clusters. 

Key words: gravitation - stellar dynamics - binaries (including multiple): close - 
stars: formation - globular clusters: general ~ open clusters and associations: general. 
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1 INTRODUCTION 

Over the last few decades, observations have revealed that 
multiple-star systems (MSSs), including binary stars and 
high-order multiples (e.g, triples, quadruples, etc.), are com- 
mon in our Galaxy. Th is is the case for the Galactic field ( e.g. 
Duguennoy fc Mayoiill991: Tokovinin 1997; Raghavan et al. 



youn g star-forming associations (e.g. iKraus et al 
2012f). old moderately-dense o p en clusters (e.g 



LathamI 120071 : iTalamantes et all l2010l : iGeller fc Mathieul 
2012f). and even ancient globular clusters (GCs) (e.g. 
Prodan fc Murrayi i201 2: Milone et al.ii2012). The most com- 



plete survey is that of lRaghavan et al.l (120101), who r e centl; 



plete survey is tnat ot IKagnavan et al.l (UUIU), wno r e cently 
updated the seminal work of iDuguennoy fc Mayoij (Il99lh 
using a volume-limited sample of field solar-type primary 
stars in the solar neighborhood selected from the Hippar- 
cos catalogue. They reported that the observed fractions of 
objects that are single, double, triple, and high-order sys- 
tems are 56 ± 2%, 33 ± 2%, 8 ± 1%, and 3 ± 1%, respec- 
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lively. Even more recently, iKraus et al.l l|2011f ) performed 
a high-resolution imaging study of the young star-forming 
region Taurus- Auriga to characterize its multiple star popu- 
lation. They found that ~ 2/3 — 3/4 of all targets are MSSs 
composed of two or more stars, while only the remaining 
~ 1/4 — 1/3 are single stars. 

The importance of dynamical encounters involving 
MSSs, even during the epoch of sta r formation, is becom- 
ing increasingly evident. iBatd l|2012f) recently presented the 
results of a statistical analysis of MSSs formed during the 
largest radiation hydrodynamical simulation of star cluster 
formation conducted to date. These models reproduce rea- 
sonably well many of the characteristics of Galactic star- 
forming regions. Interestingly, the author finds that pro- 
tostellar accretion is often terminated by dynamical in- 
teractions, and therefore that dynamics plays an impor- 
tant role in shaping the stell ar initial mass function (IMF) 
(|Podsiadlowski fc Pricelll992r ). at least in the dense clusters 
considered in this study. The importance of dynamics and 
multiplicity during the star formation process has also been 
confirmed observationally. Based on a multi-epoch search 
for wide, low-mass tertiary companions in a volume-limited 
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sample o f 118 known sgectroscopic binaries within 30 pc of 
the Sun, lAllen at alj (|2012 ) find a wide tertiary fraction of 
19.5l3y%. This is consistent with the predictions of star for- 
mation simulations, which suggest that the fraction of wide, 
low-mass companions to spectroscopic binaries is > 10%, 
and roughly twice the wide companion rate of single stars. 
This trend is thought to arise through three-body interac- 
tions during the star formation process, which transfers an- 
gular momentum away from a close pair of objects, harden- 
ing them further. 

Our theoretical understanding of how dynamical en- 
counters involving high-order multiples impact the evolu- 
tion of older star clusters is limited (e.g. iMoeckel fc Bonnelll 
[2013). The vast majority of the numerical simulations of 
star cluster evolution performed to date include only sin- 
gle stars and binaries in the initial population, and often 
these models do not allow for the dynamical formation of 
high-order multiples, thus entirely neglecting the dynamical 
input from this important stellar population. Some authors 
have recently acknowledged the need to include high-order 
multiples in their simulations, a n d list this as a future addi- 
tion ( e.g. iHvpki fc Giers3l2012l ). iGeller. Hurlev fc Mathieul 
((2013) investigate the dynamical impact of both primordial 
and dynamically formed triples on blue straggler formation 
in their Ai'-body models, and find that nearly half of the 
blue stragglers formed through collisions (between 6 and 7.5 
Gyr) resulted from stellar encounters involving hierarchi- 
cal triples. We suggest here that the inclusion of high-order 
multiples in such models may be critical for accurately re- 
producing the dynamical evolution of many star clusters as 
well as the production rates of compact binaries and stellar 
exotica. 

Recently, iLeigh fc Silla (|2011h argued that encounters 
involving triple stars should be common in the old open 
clusters NGC 188 and M67, and that this could also be the 
case for other open clusters. Importantly, however, these re- 
sults were based on the inferred properties of the binary and 
triple populations in these clusters. In this paper, we put the 
prediction of lLeigh fc Silla (|201l|) to the test using empiri- 
cally constrained binary and triple fractions in those open 
clusters for which these data are available in the literature. 
We show in Section [5] that encounters involving triples oc- 
cur roughly as often as, or even more often than, encounters 
involving single and binary stars alone in every cluster in 
our sample. More generally, using empirically-derived mul- 
tiple star catalogues for the young star-forming association 
Taurus- Auriga and the Galactic field, we show in Section [3] 
that the gravitationally-focused cross section for encounters 
increases with increasing multiplicity. This reduces the av- 
erage time between encounters and, in general, contributes 
to increasing the dynamical significance of high-order mul- 
tiple star systems. The implications of our results to various 
astrophysical phenomena are discussed in Sectional and we 
summarize the key points of this paper in Section [S] 



empirically-measured binary and triple fractions taken from 
the literature. 

iLeigh fc Silla (|201ll ) derive encounter rates for single- 
single (1-1-1), single-binary (l-|-2), binary-binary (2-1-2), 
single-triple (1+3), binary-triple (2-|-3), and triple-triple 
(3-f-3) encounters. Importantly, the ratio of any two of these 
encounter rates only depends on the binary and triple frac- 
tions, and the average geometric cross sections for singles, 
binaries and triples (e.g., the mean stellar radius, R, and 
the mean binary and tertiary semi- major axes, at, and at)- 
Therefore, the regions of parameter space where each of 
these different encounter rates dominate can easily be com- 
pared across these three different stellar populations, as we 
do in Figure [1] Each labeled polygon in this figure denotes 
the region of parameter space where the frequency of the 
given encounter type is higher than all others, and divid- 
ing lines show where the encounter frequencies are equal. 
Here we show the result using the value of at/dt for Taurus- 
Aurigqj. Perhaps the most striking result from this analysis 
is that encounters involving triples dominate over encoun- 
ters with binaries for the majority of parameter space. This 
is the case on a per star basis, even if the fraction of systems 
composed of 3 or more stars is lowQ 

In Figure [5] we show, for lines of constant binary frac- 
tion, the critical triple fraction at which encounters involving 
triples begin to dominate over encounters involving single 
and binary stars alone. To obtain each of these lines, we fix 
/i, and calculate the number of encounters involving triples 
(i.e. the sum of the numbers of H-3, 2-1-3, and 3-1-3 encoun- 
ters) over a fixed time-scale. For the same time-scale, we 
then compare this to the total number of encounters involv- 
ing binaries and singles alone (i.e. 1-1-1, 1-1-2, and 2+2). The 
critical triple fraction is shown as a function of the ratio 
dt/db- Importantly, if ft > 0.3, encounters involving triples 
will dominate over encounters involving either single or bi- 
nary stars, regardless of /;, and the ratio dt/d^. Equivalently, 
provided dt/db > 5 and ft > 0.16, triples will be undergo- 
ing encounters more often than either single or binary stars 
alone. Furthermore, for ft > 0.1, if triples are not the domi- 
nant type of interacting object, they are at least roughly as 
dynamically-active as are single and/or binary stars. 

Next we compare these analytic results to observed 
nmltiple-star populations from the literature. The points 
in Figure [l] show the o bserved binary an d triple frac- 
tions for Tauru s-Auriga IIKraus et aklboill). the Pleiade s 



(|Mermilliod et al. 1992: Bouvier. Rigaut fc Nadeau 



Praesepe (J Mermil lio d fc Mavoi 



r. nig; 
19991: 



Bouvier et al. 



1997 ). 



2001 



and the Hyades (iPatience et al.l 119981 '). We also place a 
point on Figure [l] representing a star cluster with the bi- 
nary and triple fract ions observed in the Galactic field 
(JRaghavan et al.ll20ld ). We plot these empirical results for 
illustrative purposes, but note that each survey has a unique 
level of completeness and sample selection criteria, which we 
do not attempt to correct for here. The raw numbers of sin- 
gle, binary, and triple stars are also shown for each surveyed 



2 IDENTIFYING THE DOMINANT 
ENCOUNTER TYPE 

In this section, we compare the rates of the different en- 
counter types involving singles, binaries, and triples using 



^ Specifically, at = 1773.7 AU and 05 = 130.0 AU, and assuming 
R = 1.5 R©. 

^ We define the binary and triple fractions as /;, = Nb/(Ns + Nb + 
Nt) and ft = Nt/{Ns+Nb + Nt), respectively, wliere Ns, Nb, and 
Nt are the numbers of single, binary, and triple stars, respectively 
(ignoring multiplicities higher than three in these definitions). 
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Table 1. The numbers of single, binary, and triple stars in our samples. 



Cluster Name Singles Binaries Triples 



Total 



Taurus- Auriga 

Galactic Field 

Hyadcs 

Praesepe 

Pleiades 



48 
56 
98 
43 
54 



50 
33 
59 
32 
29 



12 
8 
10 
5 
3 



110 
97 
167 
80 



^The original ^ermillio^^t_aL (1992) sample contains 88 stars, but four of these stars are associated as visual binaries, respectively. 



cluster in Table [T] Below, we briefly describe the individual 
samples. 

The most complete surveys shown here are likely those 
of the Galactic field, T aurus- Auriga, and the Hyades. The 
iRaghavan et al.l ((20101) Galactic field survey covers a com- 
plete sample of solar-type dwarfs within 25 pc of the Sun 
having 0.5 ^ [B — V) ^ 1.0 (which roughly corresponds to 
spectral types F6 - K3). This study combines multiple stars 
detected through spectroscopy, eclipsing binary surveys, sep- 
arated fringe packets (using the Center for High Angular 
Resolution Astronomy, CHARA, Array), speckle interfer- 
ometry, adaptive optics, and multiple astrometric techniques 
(e.g., astrometric orbits, common-proper motion pairs, etc.). 
This study is nearly complete for orbital sepa r ations < 10* 
AU and mass ratios > 0.1. The iKraus et al] ((201J) study 
contains all known members of Taurus- Auriga with spectral 
types between GO and M4 and masses between 2.5 M© and 
0.25 M©. Kraus et al. combine visual companions detected in 
2MASS images with companions detected through aperture 
masking with adaptive optics imagers to compile a fairly 
complete sample of multiple stars with separ ations between 
3 and 5000 AU and mass ratios > 0.1. The IPatience et al.l 
l|l998l ) Hyades speckle imaging survey covers a nearly com- 
plete sample of stars with AT < 8.5 mag (including dwarf 
stars of spectral type between about AO and K5, and some 
evolved stars), stellar separations of about 5-50 AU and mass 
ratios > 0.2. The Hyades /;, and ft values shown here also 
include multiple stars dete cted from spectrosco py and direct 
imaging (see references in lPatience et al.lll998l ). 

The Pleiades and Praesepe samples are somewhat 
less complete. For both clusters , we take the samples 
from the radial- velocity su v erys o f iMermilliod et al.l (119921 ) 
and IMermilliod fc Mayor! (|l999l ). respectively, and add 
additional compan ions detected in the adaptive o ptics 
imaging surveys of i Bouvier, Rigaut &: Nadeaul (| 19971 ) and 
iBouvier et al.l (|200ll ). respectively. For both clusters, the 
radial-velocity and adaptive optics samples do not com- 
pletely overlap, and the c overage in moderate or bital sep- 
aration is incomplete. The IMermilliod et al.l (|l992r ) Pleiades 
survey covers stars with 0.4 < (B — V) < 0.9 (about F5-K0 
spectral types), while the lBouvier. Rigaut fc Nadeaul (|l997l ) 
survey includes stars with 0.56 ^ (B — V)o ^ 1.5 (note, 
E{B — V) = 0.03) with known rotational velocities, in the 
central part of the Pleiades. The combined Pleiades radial- 
velocity and adaptive optics numbers shown in Table [T] are 
not sensitive to syst ems with separations betwee n about 5 
and 10-15 AU. The IMermilliod fc Mavorl (| 19991 ) Praesepe 
survey covers stars with 0. 4 < (B — V) < 0.8 (about F5- KO 
spectral types) , while the IBouvier et al.l (|2001 ) survey in- 
cludes proper-motion members with 0.52 < [B — V) < 1.4. 



The combined Praesepe radial-velocity and adaptive optics 
numbers shown in Table [l] are not sensitive to systems with 
separations between about 4 and 15 AU. For both clus- 
ters we also include photometric binaries identified by Mer- 
milliod et al. as being in the binary region on the color- 
magnitude diagram (brighter and redder than the single-star 
main sequence) , which may alleviate some of the incomplete- 
ness in orbital separation. 

We note that the Pleiades, Praesepe and Hyades are all 
OCs with ages on the order of several tens to hundreds of 
Myrs, and densities on the order of 1 - 10 M© pc~^. The 
ages of these clusters are comparable to the typical time 
between dynamical encounters, so that a large fraction of 
the MSSs in these clusters (and those in Taurus-Auriga) 
must be primordial. 

If these clusters have a similar ratio of at/ab as in 
Taurus- Auriga, Figure [T] shows that all of these clusters lie 
in the region where encounters involving triples will domi- 
nate over encounters involving both single and binary stars 
alone (not accounting for any incompleteness in the respec- 
tive observational surveys). 

For reference, the lGeller. Hurlev fc Mathieul (|2013l ) iV- 
body open cluster model predicts that, for dynamically 
formed solar-type triples in a rich open cluster like NGC 188, 
the ratio of at /at, will remain roughly constant over the clus- 
ter lifetime (7 Gyr for their mode l), with a mean value of 4.2 
and a standard deviation of 0.7. iGeller. Hurlev fc Mathieul 
((2013) did not include any triples in the initial population. 
Therefore their predictions are for triples formed during stel- 
lar encounters. In the Pleiades, the ratio at/at is only a 
factor of ~ 2 — 3. However the la uncertainty exceeds this 
value and, as discussed above, the Pleiades data are incom- 
plete for moderate semi-major axes. As also discussed above, 
the Taurus- Auriga data are the most complete of all of the 
open clusters in our sample, and we note that the observed 
value in Taurus- Auriga of at/ab = 13.6 is somewhat larger 
than t he theoretical prediction of lGeller. Hurlev fc Mathieul 
((2013) or the empirical value of the Pleiades. (We were un- 
able to perform a similar comparison of the at/ai, values 
for the other open clusters in our sample as these data were 
not readily available in the literature.) We conclude that the 
available data are consistent with a typical ratio at /at of at 
least a few, and possibly more. Even if at, = at, encounters 
involving triples occur roughly as frequently as encounters 
involving single stars, and encounters involving binaries oc- 
cur at most ~ 4 times as often as those involving triples. 

Thus, we conclude that encounters involving triples 
should currently be as common as encounters involving only 
single or binary stars alone in every cluster shown in Fig- 
ure [l] to within a factor of at most a few. In fact, the data 
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Figure 1. The relative rates of the different encounter types in 
the binary fraction - triple fraction plane for Taurus-Auriga, Prae- 
sepe, the Pleiades, the Hyades, and the Galactic field, as shown by 
the red points. The ratios between the geometric cross-sections for 
single, binary, and triple stars observed for Taurus-Auriga were 
used to calculate the locations of the black lines that divide the 
parameter space in the binary fraction - triple fraction plane for 
which each of the different encounter types dominates. Error bars 
on the observations denote the la uncertainties calculated using 
Poisson statistics. 



are consistent with triples being the dominant interacting 
objects in these clusters. Furthermore, given that in general 
triples are harder to detect than binaries (e.g., one may of- 
ten require both spectroscopy and direct imaging for a given 
system to detect all three components) the true locations 
of these clusters on Figure [T] may lie even further into the 
triple dominated regime. This can be properly tested when 
completeness-corrected period distributions become avail- 
able for both binaries and triples in these clusters. This will 
also allow for the results presented in this paper, which are 
calculated using average cross-sections and hence average 
encounter timescales, to be tested using timescales derived 
by integrating the full period and velocity distributions. 



3 EXTENSION TO HIGHER-ORDER 
MULTIPLICITY 

In this section, we extend our discu ssion to highe r -order mul- 
tiples using the MSS samples of iKraus et al.l l|201lh and 
iTokovininI (119971 ). These are representative of two differ- 
ent environments, namely a young (few Myrs) very low- 
density star-forming region and an older Ga l actic field pop- 
ulationlf] The MSC catalogue of iTokoviniiil iIiQQTT ) contains 



^ The field sample spans a wider range of spectral types, however 
our results remain unchanged if we limit this sample to cover the 



Figure 2. Plot showing, for lines of constant binary fraction, 
the critical triple fraction at which encounters involving triples 
begin to dominate over encounters involving singles and binaries 
alone as a function of the ratio at/ai,. Different binary fractions 
are shown using different line widths. The line width increases 
with increasing binary fraction, and the results are shown for 
ft, = 0.1,0.2,0.3,0.4,0.5. 



data on 612 (mostly) hierarchical multiple stars of multi- 
plicity three to seven, collected both from direct imaging 
and spectroscopy. The Taurus- Auriga MSS catalogue, on the 
other hand, consists of data from a high-resolution imaging 
study of 90 multiple star systems of multiplicity two to six 
with separations in the range 3-5000 AU. 

The key point we will make is that the data are consis- 
tent with a trend in which the mean cross-section for stellar 
encounters increases with increasing multiplicity. As with 
triples, this reduces the average time between encounters, 
and contributes to increasing the dynamical significance 
of higher-order multiplicity. To show this trend, we calcu- 
late for each multiplicity the gravitationally-focused cross- 
section or, more specifically, the product of the total system 
mass and the geometric cross-section. All calculations per- 
formed in this section correspond to the arithmetic mean, 
however we confirm our results using the geometric mean as 
well. All uncertainties correspond to one standard error of 
the mean|j 

First, as the number of stars increases, the geomet- 
ric cross-section (i.e. the diameter of the object) increases. 



same spectral range as the Taurus-Auriga sample. We therefore 
choose to use the complete field sample for Figure [3] 
■* We note that some of the distributions in our samples, and par- 
ticularly those of the semi-major axes, are not symmetric about 
the mean. Therefore their uncertainties about the mean are asym- 
metric. However, the sample sizes are too small to properly cal- 
culate asymmetric error bars, and therefore we simply provide 
the standard errors of the means as our best estimates of the 
uncertainties. 
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Figure 3. Plot showing the product of the mean total system 
mass (nisys; in units of Mq) and the mean maximum orbital sep- 
aration (aout ; in units of AU) , which is directly proportional to the 
mean or average gravitationally-focused cross-section for encoun- 
ters, for all types of multiples. The top inset shows these cross- 
sections for the Galactic field, for which the results are shown 
for triples, quadruples, quintuples, and sextuples. The middle in- 
set shows the cross-sections for the young star-forming association 
Taurus- Auriga, and the results are shown for binaries, triples, and 
quadruples. Finally, the bottom inset shows the cross-sections for 
the results of the hydrodynamical simulations of Bate (20f2), for 
binaries, triples, and quadruples. All error bars correspond to one 
standard error of the mean. 



and therefore the encounter time decreases. We inspect the 
MSS sample for th e star-forming association Taurus- Auriga 
IjKraus et al.ll201ll ). and compute means and standard devi- 
ations of the mean for the maximum orbital separations (i.e. 
within a given system) for all binaries, triples, and quadru- 
ples. These are (1.3±0.4) x 10^ AU, (1.8±0.5) x 10^ AU, and 
(1.4 ± 0.3) X 10^ AU, respectively. There are also two addi- 
tional multiple star systems in this sample, namely a quintu- 
ple and a sextuple, with maximum orbital separations of 3.1 
xlO^ AU, and 3.8 xlO^ AU, respecti vely. A similar a nalysis 
of the Multiple Star Catalogue (MSC: lTokovininll997l ') yields 
mean maximum (i.e. highest-order) orbital separations for 
triples, quadruples, quintuples, sextuples, and septuples of 
(l.SiO.l) X 10^ AU, (2.3±0.5) x 10^ AU, (9. Oil. 7) x 10=^ AU, 
(1.1 ±0.2) X 10* AU, and (9.4 ±3.9) x 10^ AU, respectively. 
Second, the total stellar mass tends to increase with 
increasing multiplicity, which increases the gravitationally- 
focused cross section, and thereby further decreases the 
encounter time. The mean system masses for the Taurus- 
Auriga sample (derived using th e component masses given 
in Table 5 of iKraus et al.l l201ll ) for the binaries, triples, 
quadruples, quintuple, and sextuple, respectively, are 1.0 ± 
0.1 Mq, 1.3 ± 0.1 Mq, 2.0 ± 0.3 Mq, 2.9 Mq, and 3.2 Mq. 
For comparison, in the MSC the mean system masses for 
triples, quadruples, quintuples, and sextuples, respectively. 



are4.5±0.1 Mq, 5.0±0.6 Mq, 5.6±0.6 Mq, andl3.1±1.9 
Mq (excluding systems with only minimum mass measure- 
ments) |j There is also a single septuple in the MSC sample 
with firm mass estimates for all components, yielding a total 
system mass of 28 Mq . 

The key conclusion from these results is that the 
gravitattonally-focused cross-section for encounters system- 
atically increases with increasing multiplicity. Figure [3] 
shows the product of the mean total system mass and the 
mean maximum orbital separation, which is directly pro- 
portional to the gravitationally-focused cross-section, as a 
function of multiplicity for the Galactic field (top inset) and 
Taurus- Auriga (middle inset). For comparison, we also show 
the results from the hydrodynamical simulations of iBatd 
(2012) using the values provided in their Table 3. The Spear- 
man rank correlation coefficients for these samples are 1.0 
in all cases. Thus, our proxy for the mean gravitationally- 
focused cross-section increases monotonically with increas- 
ing multiplicity in both of these observed MSS samples, 
which is als o consistent with the theoretical predictions of 
iBatd l|2012l ). 

To better quantify the statistical significance of these 
correlations, we also performed weighted lines of best-fit to 
the \og{msysaout) values sho wn in F igure [3] For the Galac- 
tic field, Taurus- Auriga, and iBatd (|2012r ) samples we find 
slopes of, respectively, 0.47 ± 0.07, 0.63 ± 0.09, and 0.40 
± 0.10. These are all inconsistent with zero at greater than 
the 3cr confidence level. To confirm this, we also calculated 
chi-squared values by comparing the log{msysaout) results 
to a fiat distribution at the respective mean for each sam- 
ple. In all three cases, we rule out the hypothesis that the 
data are consistent with a flat distribution at very high con- 
fidence (>99.99%). In summary, these data show that the 
gravitationally-focused cross-section increases with increas- 
ing multiplicity at a statistically significant level. 

The gravitationally-focused cross-section typically dom- 
inates over the geometric cross-section, especially for single 
stars and in low-velocity dispersion environments. However, 
with high-order multiples, the orbital separations can be so 
large that both cross-sections can be equally important. As 
discussed above, the data for these high-order multiples are 
also consistent with the geometric cross-section increasing 
with increased multiplicity, particularly for the MSC sam- 
ple, where the data extend from triples to a septuple. We 
note that the binding energies corresponding to the outer- 
most orbits of the MSSs in our sample are typically below 
the hard-soft boundary in their host cluster. Thus, we expect 
that most of these systems are indeed bound objects. 

The combination of the results presented in this sec- 
tion and the dynamical significance of triple stars shown in 
Section[2]suggests that high-order multiples may be very im- 
portant for many stellar dynamical processes in star clusters 
that have generally been attributed only to binary stars. We 
discuss the implications of our result below. 



■^ We note that our results remain unchanged with or without 
excluding these systems. 



© 2011 RAS, MNRAS 000, [iHl] 



6 Leigh & Geller 



4 DISCUSSION 

Multiplicity is more likely to be preserved during dynamical 
interactions in low- velocity-dispersion environments. This is 
because the hard-soft boundary is located at very long or- 
bital periods in such star clusters, so that even wide high- 
order multiples are typically classified as dynamically hard. 
We have verified this using a suite of numerical scattering 
experiments in v olving single, binary, and triple stars (see 
iLeigh fc Gellerl (|2012l ') for the details of the experiments). 
By keeping the orbital energies fixed and varying the rel- 
ative velocity at infinity, we find that decreasing the to- 
tal encounter energy on average increases the fraction of 
encounter outcomes in which high-order multiples are pro- 
duced or preserved. This higher likelihood to produce and/or 
preserve MSSs suggests that high-order multiples should 
be more common and act as more efficient heat sources in 
low-velocity-dispersion environments such as open clusters 
(OCs). 



iLeigh fc Gelleij l|2012r ) recently showed that the prob- 
ability that a direct collision will occur between any two 
stars during an encounter scales as A''^, where A'^ is the 
number of interacting stars. Therefore, on a per encounter 
basis, a higher multiplicity translates into an increased 
probability of direct stellar collisions, which may be ob- 
served as stellar exotica such as BSs. Several peculiar ex- 
amples of MSSs containing BSs have recently been identi- 
fied in OCs. Their existence is difficult to explain without 
at least some involvement from dynamical interactions. For 
example, the suspected triple sys tem S1082 in the old OC 
M67 appears to conta in two BSs (|van den Berg et al.ll200ll : 
ISandguist et al.l [20031 '). A triple system in which the outer 
co mpanion is a BS w a s rece ntly found in the OC NGC 6819 
by iTalamantes et al.l (|201G| ). A fascinati ng binary c o ntain - 
ing two BSs was also reported bv iMathieu fc Gellerl (|2009l ') 
in the old OC NGC 188. In all of these cases, mass transfer 
alone could not have produced the systems containing these 
BSs, whereas a coUisional origin or at least a subsequent 
stellar encounter(s) involving high-order multiples may po- 
tentially account for their existence. 

In fact, the multiple-star population may play a large 
role in governing both the creation and destruction rates of 
stellar exotica, in particular those th at rely on mass trans- 
fer in compact binaries. The results of lLeigh fc Gelleij (|2012l ) 
also show that physical collisions during dynamical encoun- 
ters, which become more likely with increasing multiplic- 
ity, offer a mechanism for destroying compact binaries. This 
stems from the fact that, if a collision occurs during an en- 
counter involving three or more stars, it is the most com- 
pact orbit goi ng into the interaction that w ill typically be 
destroyed fe.g. IValtonen fc Karttunenll2006l ). The details of 
the balance between creation and destruction of compact bi- 
naries resulting from encounters with MSSs in clusters is an 
important question that can be answered by the next gen- 
eration of TV-body models that i nclude significant popula- 
tions of high-order multiples (e.g. iGeller. Hurley fc Mathieul 
|2013[ ). 

On a more global scale, future A''-body models will also 
show how high-order MSSs influence the overall dynami- 
cal evolution of star clusters. For example, A'^-body sim- 
ulations have shown that the hardening of compact bina- 
ries through successive dynamical encounters in the dense 



cores of GCs can halt core collapse (e.g. lHutl ll983l). This 
evolutionary phase is often ca lled "binary burning" (e.g. 
iFregeau. Ivanova fc Rasidl2009l ). since these hardening in- 
teractions act to prevent any further increase in the central 
cluster density. However, the coUisional cross-section for a 
compact binary is very small, so that the time-scale for it 
to undergo a direct dynamical interaction with a single star 
can be very long. Consequently, the onset of the binary burn- 
ing phase is expected to e xceed a Hubble time for most o f 
the GCs in the Milky Way (JFregeau. Ivanova fc Rasidl2009l '). 
However, if most compact binaries are members of high- 
order multiples, then the time-scale for them to undergo 
direct encounters can be quite short due to the much larger 
cross-section of their parent MSS. The wide outer orbit of a 
triple effectively acts as a "net" , drawing stars in where they 
can be scattered and interact resonantly with the close inner 
binar y of the triple (JLeigh fc Sillsll201ll : iMoeckel fc Bonnelll 
120131 ). In this way, the binding energy of the compact in- 
ner binary of the triple can be tapped, and re-distributed 
throughout the cluster via relaxation processes (once the 
higher-order resonantly interacting multiple system has de- 
cayed into stable components). The role played by high- 
order multiples as heat sources in clusters should be the 
most pronounced in low-velocity-dispersion environments, 
like OCs. As discussed above, in such clusters even very 
wide orbits, which are a prerequisite for dynamical stability 
within high-order multiples, are often classified as dynami- 
cally hard. Our results suggest that triples may undergo dy- 
namical interactions roughly as frequently as binary stars, 
and hence could be relevant for the overall evolution of cer- 
tain clusters. Future detailed A''-body modeling will be nec- 
essary in order to quantify the role of triples as heat sources 
and their significance for the dynamical evolution of star 
clusters. 



5 SUMMARY 

In this paper, we address the question: How important are 
high-order MSSs to the dynamics of star clusters and the 
production/destruction of compact binaries and stellar ex- 
otica? The key conclusion is that dynamical encounters in- 
volving triple stars should be roughly as common as en- 
counters involving only single and binary stars in low- to 
moderate-density star clusters. Furthermore, the available 
data on higher-order multiples are consistent with a trend 
in which the gravitationally-focused cross-section for en- 
counters increases with increasing multiplicity. These results 
suggest that simulations of star cluster evolution should in- 
clude high-order multiples in the initial conditions and al- 
low such systems to evolve dynamically, and be created and 
destroyed, throughout the cluster evolution. This will be a 
crucial next step in the development of realistic star cluster 
simulations. 
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